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We apply the linear response theory to systems with long-range interactions described by hydro- 
| dynamic equations such as the Euler, Smoluchowski, and damped Euler equations. We analytically 

determine the response of the system submitted to a pulse and to a step function. We compare 
these results with those obtained for collisionless systems described by the Vlasov equation. We 
£SJ ' show that, in the linear regime, the evolution of a collisionless system (Vlasov) with the waterbag 

distribution is the same as the evolution of a collision-dominated gas without dissipation (Euler). 
In this analogy, the maximum velocity of the waterbag distribution plays the role of the velocity 
' of sound in the corresponding barotropic gas. When submitted to a step function, these systems 

exhibit permanent oscillations. Other distributions exhibit Landau damping and relax towards a 
\Q ■ steady state. We illustrate this behaviour with the Cauchy distribution which can be studied ana- 

lytically. We apply our results to the HMF model and obtain a generalized Curie- Weiss law for the 
magnetic susceptibility. Finally, we compare the linear response theory to the initial value problem 
' f-H '■ for the linearized Vlasov equation and report a case of algebraic damping of the initial perturbation. 

<]} PACS numbers: 

a: 
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The linear response theory is a simple and powerful tool for studying the stability of a dynamical system and 
its response to an external perturbation. Kubo [l| used it to determine the variation of the average value of an 
"observable" due to an applied "force" . The linear response theory found a lot of applications in statistical mechanics 
and kinetic theory , and in the theory of simple liquids Q . Its most important applications concern the derivation 
' O I of expressions for the transport coefficients of hydrodynamics, including the electrical and thermal conductivity, and 
. the mobility of a Brownian particle immersed in a fluid. 

The linear response theory has also been applied to long-range interacting systems 043- In particular, it has 
been used to compute the friction force experienced by a star passing through a globular cluster [8l-fl(j| or the drift 
of a point vortex moving in a background shear [TlL 12 1. These results can be generalized to other systems with 
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I. INTRODUCTION 



long-range interactions [13|. Recently, the linear response theory has been applied to the Hamiltonian Mean Field 
(HMF) model [l4[ in order to determine how a long-range system in a quasistationary state (QSS) responds to an 
external perturbation [l5l [l6j . In the examples discussed above, the linear response theory is based on the Liouvillc 

■ equation for a Hamiltonian system of particles in interaction, or on the Vlasov equation for a collisionless gas. In the 
I/"") ' present paper, we apply the linear response theory to a collision-dominated gas described by hydrodynamic equations 

such as the Euler, Smoluchowski, and damped Euler equations. Specifically, we consider the response of this system 

■ to a weak external perturbation, and we compare the results with those obtained for a collisionless gas. 

£SJ , The paper is organized as follows. In Section [Hi we recall general relations satisfied by the polarization and 
7-H ■ response functions. As an application, we consider the response of the system to a pulse and to a step function. In 
Section IIII1 we apply the linear response theory to a collision-dominated gas described by hydrodynamic equations. 
• i-h ■ For the sake of generality, we consider the damped Euler equation including a friction force — £u proportional to 
1 the velocity. For £ = 0, we recover the pure Euler equation describing an ideal gas, and for £ — > +oo, we get 
^ . the Smoluchowski equation describing the overdamped motion of Brownian particles. Explicit expressions of the 
response and polarization functions are obtained for these systems. In Section IIV1 we compare these results with 



those obtained for a collisionless gas described by the Vlasov equation. We show that, in the linear regime, the 
evolution of a collisionless system (Vlasov) with the waterbag distribution is the same as the evolution of a collision- 
dominated gas without dissipation (Euler). In this analogy, the maximum velocity of the waterbag distribution plays 
the role of the velocity of sound in the corresponding barotropic gas. When submitted to a step function, these systems 
exhibit permanent oscillations. In Section [Vj we treat the case of the Cauchy distribution for which analytical results 
can be obtained. When submitted to a step function, this distribution exhibits Landau damping and relaxes towards 
a steady state. In Section I VII we study the evolution of the perturbed density in the linear regime and establish 
hydrodynamic equations for a collision-dominated gas and a collisionless system. In Section IVII1 we consider the 
asymptotic distribution of a stable system submitted to a step function and show that it coincides with the steady 
distribution of this system under a weak external field. In Section IVIII1 we apply our general results to the HMF 
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model and consider specific distribution functions such as the isothermal, polytropic, waterbag, and Fermi-Dirac 
distributions. We show that the magnetic susceptibility of the equilibrium state is given by a generalized Curie- Weiss 
law. In Section IIX1 we compare the response of a collisionlcss system submitted to a weak external potential with 
the evolution of an isolated system submitted to an initial disturbance (the so-called initial value problem of the 
linearized Vlasov equation) and we report a case of algebraic damping. In Appendix we discuss the connection 
between a collisionless system described by a distribution function depending only on the individual energy and the 
corresponding barotropic gas. In Appendix [B] wc illustrate some results for Coulombian plasmas, self-gravitating 
systems, and for the HMF model. 

II. POLARIZATION AND RESPONSE FUNCTIONS 
A. General results 

We examine the response of a spatially homogeneous system at equilibrium to a small external potential $ e (r, f). 
We follow the presentation given by Binney and Tremaine . The perturbation caused by the external potential may 
be written as 

<y$(r,t)= f dr'dt' R(r-r',t-t')$ e (r',t'), (1) 

where i?(x, r) is the response function. Causality requires that i?(x, r) = for r < 0. We introduce the polarization 
function P(x, t) which relates the perturbation to the total potential: 

., W h/« P (->flW,7, t » W )l. (2) 

Once again, P(x, r) = for r < 0. The polarization function corresponds to the response of the system when collective 
effects arc neglected. Indeed, if we neglect the self-interaction of the system 8$ in the calculation of the response 
function R (see below), we obtain the polarization function P. 

Since the external potential is introduced at t = (say), it is convenient to use Laplace transforms in time and 
Fourier transforms in space. The Fourier-Laplace transform of the perturbed potential <5$(r,t) is defined by 

/dr f + °° 
j^S I dte-^-^S^t). (3) 

This expression for the Laplace transform is valid for Im(cj) sufficiently large. For the remaining part of the complex 
lo plane, it is defined by an analytic continuation. The inverse Fourier-Laplace transform is 

5$(r,t)= J dk J ^e i ( k r -^ ( 5l'(k,w), (4) 

where the Laplace contour C in the complex oj plane must pass above all poles of the integrand. Taking the Fourier 
transforms of Eqs. ([1]) and (|2|), wc obtain 

8&(k,t) = / dt' R(k,t-t')3> e (k,t'), (5) 



5$(k,t)= [ dt'P(k,t-t') [$ e (k,t') + £$(k,t')l . (6) 
Jo L - 1 

Taking the Fourier-Laplace transform of Eqs. ([1]) and ([2]), we get 

a$(k,«) = .R(k J w)*e(k > w) J (7) 

(5$(k, lu) = P{k, uj) [$ e (k, uj) + S$(k, uj)} . (8) 
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The Fourier-Laplace transforms of the response and polarization functions are related to each other by 

, P(k,u) „ . RCk,u) 



If we define the dielectric function by 



we obtain 



c(k,«;) = l-P(k > w) > (10) 



un x l-e(k,u;) 

R(Ku) = — jt- — — . 11) 
e(k, uj) 

B. The response to a pulse 

We consider the response of the system to a "pulse" localized at t = 0. It can be represented by a Dirac distribution 

$ e (M)=(f(t). (12) 

We have directly written the Fourier transform of the external potential and, for simplicity, we have assumed that 
the perturbation is independent on the mode k (this situation can be straightforwardly generalized by multiplying 
the results by some amplitude <J? e (k)). The Laplace transform of the external potential is 

* e (k,w) = l. (13) 

According to Eqs. (JT]) and ([IB"]) , the perturbation caused by a pulse is equal to the response function: 

5$(k,w) = R(k,ui). (14) 

Taking the inverse Laplace transform of this expression, and using Eq. (|11[) . we obtain 

,4(M) = / (15) 

J c 2tt e(k,w) 

The poles of the integrand correspond to the complex pulsations ui a (k) for which the dielectric function vanishes: 
e(k, cj a (k)) = 0. This defines the dispersion relation. The evolution of the perturbation depends on the position of 
the zeros of the dielectric function in the complex plane. Using the Cauchy residue theorem, we have 



e -iu a (k)t 



Res 1 -^ 



e (k,w) 



(16) 

l>— cj a (k) 



where the sum runs over the whole set of poles and we have assumed, for simplicity, that the singularities are simple 
poles. In the following, we shall omit the subscript a for brevity. If at least one zero tu of the dielectric function lies 
in the upper half plane (i.e. uJi > 0), the system is unstable, and the perturbation grows exponentially rapidly with 
the rate (u>i) max corresponding to the zero with the largest value of the imaginary pulsation. If all the zeros u> of the 
dielectric function strictly lie in the lower half-plane (i.e. UJi < 0), the system is stable, and the perturbation decays 
to zero exponentially rapidly with the rate |wj| m in corresponding to the zero with the smallest value of the imaginary 
pulsation in absolute value. Finally, if the zeros lie on the real axis (i.e. w, =0), the system is marginally stable and 
the perturbation displays an oscillating behavior with the pulsation cj r . For more details, we refer to 0, H|. 

C. The response to a step function 

We consider the response of the system to a constant potential that is "switched on" suddenly at t = 0. It can be 
represented by a step function 

$e(k,i) = F(t)$ e (k), (17) 
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where H(t) — for t < and H(t) = 1 for t > (Heavisidc function). The Laplace transform of the external potential 



is 



* e (k,w) = -$ e (k). 



(18) 



According to Eqs. (J7J and (jlijp. the perturbation caused by a step function is 



<J$(k,w) - iE(k,w)-* e (k). 



(19) 



Taking the inverse Laplace transform of this expression, and using Eq. ([lip , we obtain 



1 — e(k, w) i 
e(k, w) w 



(20) 



The integrand presents a pole at w = that gives rise to a constant term (1 — e(k, 0))/e(k, 0). On the other hand, 
the temporal evolution of the perturbation depends on the position of the zeros of the dielectric function e(k, u>) in 
the complex plane. Using the Cauchy residue theorem, we have 



<y$(k,t) l-e(k,0) 



*e(k) 



e(k,0) 



E 1 



Res 



l-e(k,cj) 
e(k, u>)uj 



(21) 



where we have assumed, for simplicity, that the singularities are simple poles. If at least one zero u> of the dielectric 
function lies in the upper half plane (i.e. Wj > 0), the system is unstable, and the perturbation grows exponentially 
rapidly with the rate (wi) max . If all the zeros u> of the dielectric function e(k, u>) strictly lie in the lower half-plane 
(i.e. u>i < 0), the late time evolution of the perturbation <5$(k, £) is dominated by the pole at ui — 0. As a result, the 
perturbation tends to the asymptotic value 



{ *~< k > - i i(^f *•*>• 



(22) 



for t — > +oo. Finally, if the zeros lie on the real axis (i.e. = 0), the perturbation displays an oscillating behavior 
about the value (|22[) with the pulsation u> r . 



III. THE DAMPED BAROTROPIC EULER EQUATION 



We consider a system of particles with long-range interactions interacting via a binary potential u(\r — r'|). We 
assume that the particles also experience collisions due to short-range interactions. In the absence of dissipation, this 
system is described by the Euler equations. These hydrodynamic equations are the correct description of a collision- 
dominated gas. For the sake of generality, we also allow for the possibility that the particles of the system move in an 
inert medium and experience a friction force proportional to their velocity. This is the case, for example, for colloidal 
particles immersed in a fluid. In order to study these two situations in a unified manner, we consider the damped 
Euler equations fi~7j : 



dp 
dt 



V ■ (pu) = 0, 



(23) 



du 

~dt 



+ (u • V)u 



-Vp - pV$ - pV<I> e - £pu, 



(24) 



$(r,i) = J u(\r-r'\)p(r',t)dr', (25) 

where £ is the friction coefficient. We have distinguished the potential $(r, f) produced by the particles from the 
external potential $ e (r, i). To close the system of equations, we need to specify the equation of state. We consider a 
barotropic gas in which the pressure is a function of the density: p(r, t) = p[p{v, t)]. The specification of the equation 
of state p = p(p) completely closes the system of equations. 
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The friction coefficient £ measures the strength of dissipative effects. For £ = 0, Eqs. (|23|) - (l25l) reduce to the Eulcr 
equations for a perfect fluid. Alternatively, in the strong friction limit £ — > +00, we can neglect inertial effects in Eq. 
pi]) and obtain 

pu~ -i(Vp + /9V$ + pV$ e ). (26) 

Substituting this relation in the equation of continuity (|23|) , we get the generalized Smoluchowski equation [l7j : 

dp 



6t V 



(27) 



$(r,i) = / u{\r-r'\)p(r',t)dr'. (28) 



This equation describe the dynamics of Langevin particles in interaction in a strong friction limit. Therefore, the 
damped Eulcr equation makes the link between the Euler equation (£ = 0) and the Smoluchowski equation (£ — > +00). 
For a more detailed discussion of these equations, we refer to fl8j . 

When the potential of interaction u(\r — r'|) is the gravitational potential, the mean field equation (f2"5j) or (j2"8"]) 
reduces to the Poisson equation. The Euler-Poisson system describes a self-gravitating collision-dominated gas like a 
barotropic star for example [|| . The damped Euler-Poisson system may describe the dynamics of dust particles in the 
solar nebula [l9| . In the strong friction limit, the generalized Smoluchowski-Poisson system describes self-gravitating 
Langevin par ticles [20j . Finally, the case of a cosine potential of interaction in ID has been considered in relation to 
the HMF [Hi and BMF (p. 86 of (HJ) models. 

A. The linearized damped Euler equations 

We consider a system in a steady state with uniform density p and examine its response to a small external 
potential $ e (r, i). Since the perturbation is small, we can develop a linear response theory. The linearized damped 
Eulcr equations are 

^+pV-u = 0, (29) 

9 ~E = ^ v5p ~ pv6q> ~ pv ® e ~ (30) 

5$(r,t) = J u(\r-r'\)8p(r',t)dr', (31) 

where c 2 s = p'(p) is the velocity of sound in the gaseous medium. Taking the Fourier-Laplace transform of these 
equations, we get 

— iLuSp + ipk ■ u = 0, (32) 

— iujph = —Cgi5pk — ipS&k — ip$ e k — £pu, (33) 

5$ = (2n) d u(k)Sp. (34) 
Taking the scalar product of Eq. ([33| with k, and using Eqs. ([32]) and (|34| . we find that 

<5l>(k,w) = i?(fc,w)$ e (k,a;), (35) 

with the response function 

Mk \ = {2*) d u{k)e P 
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Using Eqs. (j9)) and (jTUJ) , we obtain the polarization function 



w z + iciu; — c^k 



and the dielectric function 



+ lt;UJ — CgK z 



B. The solution of the dispersion relation 

The dispersion relation e(k, uS) = can be written explicitly as 

u> 2 + i& - u>%(k) = 0, (39) 

where uj 2 (k) = c 2 s k 2 + (2tt) q l u{k)k 2 p. The dispersion relation determines the complex pulsation u> as a function of the 
wavenumber k. We find that 



u _ (40) 

where A(fc) = 4wQ(fc) — £ 2 . If cJq(&) > 0, the imaginary part of the complex pulsation is negative, implying stability. 
If Wg(fc) < 0, the imaginary part of the complex pulsation with the sign + is positive, implying instability. Therefore, 
the system is stable with respect to a perturbation with wavenumber k when 

e (M) = i + ^^>o > (41) 

and unstable otherwise. This stability criterion can also be obtained from the Nyquist theorem [22j |. It provides a 
generalization of the Jeans instability criterion in astrophysics @. For repulsive potentials u(k) > (as in plasma 
physics), the system is always stable. For attractive potentials u(k) < (as in astrophysics), the system may be 
unstable to some wavelengths. 

A mode is a perturbation that can be sustained without external forces. It is therefore the solution of Eqs. (f29)) - 
(l3"TT) with $ e = 0. The modes are usually of the form Sp oc e *(k-p— «(fc)t) w here us(k) is the solution of the dispersion 
relation (p9|) . When A(fc) 7^ 0, 5p is a linear combination of the two modes w±(fc) given by Eq. (|40|) . When 
A(fc) = 0, w+(k) = uj-(k) = and the solution of Eqs. (f29|) -(|3T |) is a linear combination of Sp oc e ~£*/ 2 e lk ' r 

and Sp oc te~^*/ 2 e 4k ' r (this may be seen, for example, on the hydrodynamic equation (|113p that is equivalent to Eqs. 

©-(Si)). 

For the Euler equation (£ = 0), the dispersion relation reduces to lu 2 = w 2 (fc). In the stable case, the perturbation 
oscillates with a pulsation ui r = ±oj (fc)- In the unstable case, the perturbation grows exponentially rapidly with 
a growth rate uii = y— Wg(fc) > (the second mode is damped exponentially rapidly with a damping rate uii = 
-V^p)<0). 

For the Smoluchowski equation (£ — > +00), the dispersion relation reduces to i£cj = Wg(/c). In the stable case, the 
perturbation is damped exponentially rapidly with a damping rate cjj = — Wg(fc)/^ < 0. In the unstable case, the 
perturbation grows exponentially rapidly with a growth rate = — cjQ(fc)/£ > 0. 

We now consider the damped Euler equation. In the stable case, we have to distinguish three cases: If A(fc) > 0, 
the perturbation oscillates with a pulsation uj r — A(fc)/2 while being damped exponentially rapidly with a 
damping rate = — £/2 < 0; if A(fc) < 0, the perturbation is damped exponentially rapidly with a damping rate 
uji = (— £ + y— A(fc))/2 < (the second mode is damped more rapidly at a rate = (— ^ — y/—A(k))/2 < 0); if 
A(fc) = 0, the temporal evolution of the perturbation behaves as te^^l 2 . In the unstable case, the perturbation grows 
exponentially rapidly with a growth rate uii = (— £ + \J — A(fc))/2 > (the second mode is damped exponentially 
rapidly with a damping rate = (— £ — y / — A(fc))/2 < 0). 

The dispersion relation (|39|) has been studied in detail for specific potentials of interaction (self-gravitating systems, 
chemotaxis, plasmas, and BMF model) in [23l |24j. 
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C. The response to a pulse 



The response of the system to a pulse is given by Eq. (fl5j) . Using Eq. (j38[) , the evolution of the perturbation can 
be written as 

5$(k, t) = {2ir) d u(k)k 2 pl(k, t), (42) 

with 

I(k,t)= [ IV^-j— i jr^. (43) 

We recall that these equations determine the response function R(k,t) [see Eq. (fl"4)) ]. The integral (|43)) can be easily 
calculated with the residue theorem. 



1. Stable case 

We first consider the stable case Wo(k) > 0. For the Euler equation, 



For the Smoluchowski equation, 



For the damped Euler equation, 



/(M) = -^#. (44) 



1 "0< fc > t 

I{k,t)=--e — *. (45) 



I(k,t) = -^= e -^sin Qv/AW*) , (A(*) > 0), 



(46) 



: e -«*/ 2 sinh Q,/=A(fc)t) , (A(fc) < 0), (47) 



2 



I(k,t) = -te' it/2 , (A(fc) = 0). (48) 
The evolution of the perturbation is consistent with the discussion given in Section IIIIBI 



2. Unstable case 



We now consider the unstable case uj 2 (k) < and define 7o(fc) = —u> 2 (k). We can either compute the integral (j43 
with the residue theorem or replace too(k) by ijo(k) in Eqs. (|4"4"]) -(|48 |) . For the Euler equation, 



/(M) = _smh[^] 
7o (fc) 



For the Smoluchowski equation, 



1 ^i') , 

I(k,t) = --e- a t- t . (50) 



For the damped Euler equation, 

2 



I(M) = - 



-^sinhQv^ATfc)^ , (51) 



with A(fc) = — £ 2 — 47 2 (fc) < 0. The evolution of the perturbation is consistent with the discussion given in Section 

mm 
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D. The response to a step function 



The response of the system to a step function is given by Eq. ([20]) . Using Eq. (|38j) . the evolution of the perturbation 
can be written as 



<5$(k,i) = (2ir) d u(k)k 2 P J(k,t)$ e (k), 



with 



J(k,t) 



du> 



1 



i c 2ir cu 2 + i£uj — 0J 2 (k) ui 
The integral (j5"3")l can be easily calculated with the residue theorem. 



(52) 



(53) 



1. Stable case 



We first consider the stable case oj 2 {k) > 0. For the Euler equation, 



For the Smoluchowski equation, 

For the damped Euler equation, 
J(k,t) 



J(k,t) 



J(k,t) 



u 2 (k) uj 2 (k) 



-it/2 



j 2 (k) 



1 - cos[u (k)i\). 



1 — e I 



VW) V 2 / V 2 > 



(A(k) > 0), 



(54) 



(55) 



(56) 



J(k,t) 



-it/2 



UJ 2 (k) uj 2 (k) 



v/=A(fe) 



: sinh 



cosh 



(A(fc) < 0), (57) 



J(M) 



4 

e 



-it/2 



e 



(A(fc) = 0). 



(58) 



The pole at u> — contributes to the integral J(k,t) by a constant term l/<Jo(fc). When £ = 0, the zeros of 
the dielectric function e(k,ui) lie on the real axis (w r = ±u>o(k), tOi = 0). As a result, the perturbation oscillates 
indefinitely about l/uj 2 (k) with a pulsation ±u>o(k) [see Eq. (l54l) ]. When £ > 0, the zeros of the dielectric function 
e(k,u) strictly lie on the lower half plane (cjj < 0). In that case, the perturbation asymptotically tends to l/oj 2 (k) 
for t -> +oo [see Eqs. ([55} -([58 ]l ]. 



2. Unstable case 



We now consider the unstable case uJ 2 (k) < and define 7q(/c) = —u> 2 (k). We can either compute the integral (f5"3")) 
with the residue theorem or replace wo(fc) by ijo(k) in Eqs. (l54l) - ([58|) . For the Euler equation, 



For the Smoluchowski equation, 



J(k,t) 



J{k,t) 



7o 2 (fc) 



(1 -cosh[ 7o (fc)t]). 



7o 2 (fc) 



1 - e « 



(59) 



(60) 
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For the damped Euler equation, 

1 .... r ? i j~\(h\ \ ( j-\(k\ x 

(61) 



J(k, t) — 2 2 e ^*/ 2 

7o0) 7o( fc ) 



sinh ( VEM t + C0S h [ 



v^A(fc) 



with A(/c) = -i 2 - 4 7 2 (fc) < 0. 

Remark: In the unstable case, the perturbation grows exponentially rapidly as explained in Section flll Bl Of course, 
the linear response theory ceases to be valid when the perturbation has grown significatively, so the expressions 
obtained in Sections MI C 21 and MI D 21 are only valid for sufficiently short times. 

E. The polarization function 

The polarization function is given by Eq. (|37[) . It can be written as 

P(k,t) = {2Tr) d u{k)k 2 pK{k,t), (62) 

with 

K(k,t) = f 2T2- ( 63 ) 

Comparing the polarization function (|37p with the response function Q36p , we see that they only differ by the replace- 
ment of u>o(k) by c 2 s k 2 (this amounts to neglecting the self-interaction as indicated in Section [TTJ) . We can therefore 
readily adapt the results of Section MI CI For the Euler equation, 

c s k 



For the Smoluchowski equation, 



For the damped Euler equation, 



K(k,t) = ~e-^f t . (65) 



K(k,t) = -_^= e -««/a sin Q^/a^)*), (A(fc) > 0), (66) 
^__, e -«*/ 2 sinh Qv^A(fc)t) , (A(fc) < 0), 



K(k,t) = -^—-e- (t/2 sinh ( ^^A(fc)t ) , (A(fc) < 0), (67) 



K{k,t) = -te^ t/2 , (A(ifc) = 0), (68) 

with A(k) = 4c 2 k 2 - C 2 . 

IV. THE VLASOV EQUATION 

We consider a system of particles with long-range interactions interacting via a binary potential u(\r — r'|). We 
assume that "collisions" (correlations, grainincss, finite N effects) are negligible. In that case, the system is described 
by the Vlasov equation 

| + v.|-(v $ + va> e ).^ = o, (69) 

*(r,t)= / u(\v-v'\)p(r',t)dv', (70) 
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where /(r, v, t) is the distribution function and p(r,t) = J /(r, v, t)dv the density. As before, we have distinguished 
the potential $(r,£) produced by the particles from the external potential $ e (r,i). The Vlasov equation, which is 
based on a mean field approximation, describes the collisionless evolution of stellar systems, plasmas, and of the HMF 
model. It is rigorously valid for systems with long-range interaction in a proper thermodynamic limit N — > +00 j25j . 
Systems with long-range interactions are known to organize spontaneously into quasi stationary states (QSSs) that 
are steady states of the Vlasov equation @, 0] ■ Galaxies in astrophysics and large-scale vortices in 2D hydrodynamics 
are examples of such QSSs [26[. These QSSs have also been studied extensively for toy models like the HMF model 
[l4l [27| - |46| or simplified models of gravitational dynamics [47l-f59j . The linear response theory may be applied to 
systems in such QSSs. For simplicity, we restrict ourselves to spatially homogeneous systems as in [15| . The case of 
spatially inhomogeneous systems can be treated with angle-action variables as in (l6j . 

A. The linearized Vlasov equation 

We consider a system in a steady state with uniform density p and distribution function /(v), and examine its 
response to a small external potential <j> e (r,t). Since the perturbation is small, we can develop a linear response 
theory. The linearized Vlasov equation is 

g + T .g_ ( v ( . + v..,.f£-o, (n, 

5§(r,t)= f u(\r-r'\)5p(r',t)dr'. (72) 
Taking the Fourier-Laplace transform of this equation, and assuming Sf(t = 0) = 0, we get 



k 



Of 



Sf(k,v,cj) = - <5$(k, w ) + $ e (k, w ) , (73) 

K ' V CJ L 

<5$(k,w) = (2ir) d u(k)5p(k,uj). (74) 

Integrating Eq. ([73f over the velocity, and using Eq. (|74| , we find that the Fourier-Laplace transform of the perturbed 
potential and of the perturbed distribution function are given by 

6*Qt,«j)= 1 ~' (k ' Cj) $ e (k,a;), (75) 
e(K, Lo) 

k . &£ 1 

<J/(k, v, u) = *!_- __$ e (k, u), (76) 

k • v — uj e(k, uj) 

with the dielectric function 

f k- ^ 

e(k, u) = 1 - (2ir) d u(k) / dv. (77) 

We consider a c?-dimensional space. It is convenient to take the Ud-axis in the direction of k. If we integrate over 
V1...V4 and note v for Vd and f(v) for J f dvi...dvd-i (in the following, f(v) will be called the reduced distribution 
function), we obtain 

e(k,oj) = l-(2ir) d u(k) [ -J-^—dv, (78) 
J L v- uj/k 

where the integration has to be performed along the Landau contour L [60| . For the Maxwell distribution (| 145[) , we 
can write the dielectric function in the form 

e(k, w) = 1 + (2n) d u(k)(3pW (v^|) , (79) 



where 

— 

v27T J L x z 

is the plasma dispersion function [6l| . 



W(z) = -£= I ^—e- x2 ' 2 d Xl (80) 
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B. The dispersion relation 

The dispersion relation e(k, uS) = can be written explicitly as 



\-{2^) d u{k) [ f ' i - V \ l dv = 0. (81) 
v — u/k 



L 



It is in general difficult to solve this equation analytically except for the waterbag and the Cauchy distributions (see 
Sections I IV El and [V} . Analytical results can be obtained for the Maxwellian distribution in some asymptotic limits 
0;SH3- I n general, the imaginary part of the complex pulsation ui is non-zero, leading to Landau damping (w.j < 0) 
or Landau growth (w, > 0). Using the Nyquist theorem [22|], it is possible to obtain a general criterion of dynamical 
stability. If the distribution function f(v) has a single maximum at v = 0, it can be shown that the system is stable 
with respect to a perturbation with wavenumber k when 

r + °° f(v) 

e(fc, 0) = 1 - {2v) d u{k) / LU. dv > 0, (82) 

J-oo v 

and unstable otherwise. This is a generalization of the Jeans stability criterion in astrophysics @. For repulsive 
potentials u{k) > (as in plasma physics), a single humped distribution is always stable. For attractive potentials 
it(k) < (as in astrophysics), a single humped distribution may be unstable to some wavelengths. When the 
distribution function is of the form / = F(v 2 /2), using the notion of "corresponding barotropic gas" (see Appendix 
|A"|) , one can show that the stability criterion ([52")) can be written as 

^0)^l+ (27r) ? fcV >0, (83) 



where c s is the velocity of sound in the corresponding barotropic gas. The equivalence between Eqs. (|82p and 
is due to the identity (|A3[) . On the other hand, Eq. ([55)) is equivalent to the stability condition (|4"T)) for a collision- 
dominated gas. Therefore, for spatially homogeneous distributions, a collisionlcss s yste m with a distribution function 
/ = F(v 2 /2) is stable if, and only, if the corresponding barotropic gas is stable [22|]. This equivalence is not true 
anymore for spatially inhomogeneous distributions 0, [fi2j . 

C. The polarization function 

The Fourier-Laplace transform of the polarization function is related to the dielectric function by Eq. (|10|) . Using 
Eq. ([771), we obtain 

r k- 2£ 

P(k,w) = (27r) d w(fc) / ^dv. (84) 

J k • v — uj 

Taking the inverse Laplace transform of this equation, we get 

P(k, t) = (27t) d u(k)i J k • |^e- lk ' v * dv. (85) 
For the Maxwell distribution (|145[) . we have the explicit result 

P(M) = -(2n) d u(k)k 2 pte-*W-. (86) 

The polarization function of a collisionless system with a Maxwell distribution decays to zero, contrary to the polar- 
ization of an ideal fluid which oscillates [see Eq. ([Ml) ]. This is a consequence of phase mixing @. 

D. An integral equation for the response function 

The Fourier transform of the linearized Vlasov equation (|71[) is 

^ + ik ■ vSf -ik~ [<J*(k, f) + & e (k, t)] = 0. (87) 
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Assuming Sf(t = 0) = 0, Eq. (|87|) can be integrated into 



5f(k,v,t)=ik-^- dt' 



5${k,t') + $ e (k,t') 



-ik-v(t-t') 



where we recall that <5$(k, t) depends implicitly on 6f(k, v, t) through Eq. (I7H) . Integrating Eq. 
we get 



Softs., t) = i I dv k ■ ^ / dt' 
ov Jo 



5$(k,i') + $e(M') 



-ik-v(t-t') 



Using Eq. (j74p , we find that the perturbed potential satisfies an equation of the form 
6$(k,t) = (2Tr) d u(k)i f dvk~ [ dt' \s$(k,t') + $ e (k,t' 



-ik-v(t-t') 



over the velocity, 
(89) 

(90) 



Comparing Eq. (|90[) with Eq. ([5]), we recover the expression (|85[) of the polarization function. On the other hand, 
for a pulse $ e (k, t') = S(t'), the perturbed potential <5$(k, t) is equal to the response function R(k,t) [see Eq. (fT4]) ]. 
Therefore, the response function satisfies the integral equation 



R(k,t) = (27r) d u(k)i J dvk 



dv 



(2 



n) d u(k)i Jdvk-^-J^ d^'e- <k • v ( t - t ' ) i^(k,^ , ). 



For the Maxwellian distribution, we can reduce the foregoing equation to the form 



R(k,t) = -(2ir) d u(k)k 2 pte-^ - (2i:) d u(k)k 2 p I dt' (t - t^e"^^- 4 ') 2 R(k,t') 



(91) 



(92) 



The first term accounts for phase mixing [see Eq. (1861) ] and the second term for Landau damping or Landau growth 



E. The waterbag distribution 

We consider a reduced distribution function of the form f(v) = p/2v m for — v m < v < v m and f(v) = otherwise 
(we note that f(v) — > pS(v) when v m —> 0). This is the so-called waterbag distribution. Using f'(v) = (p/2v m )[6(v + 
v m ) — S(v — Vm)], the dielectric function (|78j) is explicitly given by 



(2n) d u(k)Pp 
s(k,u) = 1-— 5 "TTT2 • 



CJ 2 — w^fc 2 

The dispersion relation e(k,ui) = can be written as 

u? = v m k 2 + (2ir) d u(k)k 2 p. 



(93) 



(94) 



The pulsation is purely real or purely imaginary. The system is stable with respect to a perturbation with wavenumber 
k when 



„ . (2ir) d u(k)p 



(95) 



and unstable otherwise. In the stable case, the perturbation oscillates with a pulsation uj r = ±\/ljj 2 without being 
damped. In the unstable case, the perturbation grows exponentially rapidly with a rate = \/—oj 2 > (the second 
mode is damped exponentially rapidly with a rate cj; = — y— w 2 < 0). 

The important point to notice is that the dielectric function (|93| of a collisionless system described by the Vlasov 
equation with the waterbag distribution coincides with the dielectric function (|38p of a collision-dominated gas de- 
scribed by the Euler equation (£ = 0), provided that the velocity of sound c s is replaced by the maximum velocity v m . 
Therefore, the results obtained in Section HLT1 for the Euler equation are immediately applicable to the Vlasov equation 
when the unperturbed system is described by the waterbag distribution. A property of the waterbag distribution is 
that it does not experience Landau damping so the perturbation has a purely oscillatory behaviour in the stable case 
as for a perfect gas [see Eqs. (|4"4")l . (|54p . and (|64p] . 
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THE CAUCHY DISTRIBUTION 



It is in general difficult to solve the dispersion relation (|8Tj) of the linearized Vlasov equation analytically, even for 
the Maxwell distribution. Analytical results can be obtained for the waterbag distribution (sec Section IIVE[) , but 
this distribution is very particular because the zeros of the dielectric function lie on the real axis (in the stable case) 
so the perturbations do not experience Landau damping. This is to be contrasted with most distribution functions, 
including the Maxwellian. In this respect, it can be interesting to consider the Cauchy distribution 

/(«) = -£-— ^ (96) 
™° 1 + k 

which is less peculiar than the waterbag distribution, and for which analytical results can be obtained (a drawback 
of this distribution is that the mean square velocity diverges). 

A. Dispersion relation 

For the Cauchy distribution, the dielectric function ([75]) can be written as 



2p r x 



e(k,uj) = 1 + (2n) d u(k)-^ / 7T —j^ nz — dx. (97) 



For uj lying in the upper half plane, the integral can be computed analytically by adding to the real axis a large 
semi-circle in the lower half-plane and using the Cauchy residue theorem for a function with a double pole at x = —i. 
The evaluation of the integral leads to the result 

. (2it) d u(k)pk 2 , . 

e(k,u>) = l- \.J ( '\ 2 . (98) 

This result is then extended to any ui by analytic continuation. The dispersion relation e(k,uj) = is readily solved, 
giving the two complex roots: 

uj = ±^(2n) d u{k)k 2 p - iu k. (99) 

For a repulsive potential u{k) > 0, the system is stable. The perturbation oscillates with a pulsation ui r = 
±y (2n) d u{k)k 2 p and is damped exponentially rapidly with a rate = — uok < 0. For an attractive potential 
u(k) < 0, the system is stable with respect to a perturbation with wavenumber k when 

„ , (2n) d u(k)p , . 

e(k, 0=1 + ^ — ' 9 V IH > 0, (100 

«o 

and uns table otherwis e. In the stable case, the perturbation is damped exponentially rapidly with a rate 
u>i = y {2ir) d \u(k)\k 2 p — u$k < 0. In the unstable case, the perturbation grows exponentially rapidly with a 
rate cjj = y {2ir) d \u(k)\k 2 p — uok > 0. The other mode is always damped exponentially rapidly with a rate 
= -^/(2Tr) d \u(k)\k 2 p-u k < 0. 

We note that the dielectric function (|9"8"|) of a collisionless system described by the Vlasov equation with the Cauchy 
distribution is similar to the dielectric function (|38|) of a collision-dominated gas described by the damped Eulcr 
equations (£ > 0), where Uq plays the role of c s and 2fcu the role of £. They are, however, different because the 
friction coefficient £ in Eq. (|38|) is constant while the damping term 2kuo in Eq. (|98|) depends on the wavenumber k 
(they coincide only for the HMF model for which k = 1). 

B. Polarization function 

The Fourier-Laplace transform of the polarization function is related to the dielectric function by Eq. (|10[) . Using 
Eq. ([98|), we get 

(2n) d u(k)pk 2 nnn 
P{k,w) = -f7i : — S2~- (101) 

[ikUQ + U!) Z 
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Taking the inverse Laplace transform of Eq. (|101|) , we obtain the explicit result 

P(M) = -(2ir) d u(k)pk 2 te- Uokt . 



(102) 



The polarization function decreases exponentially rapidly as a result of phase mixing (we note that the polarization 
function associated with the Maxwellian distribution has a Gaussian decay). 



C. The response to a pulse 

The Fourier-Laplace transform of the response function is related to the dielectric function by Eq. (jlll) . Using Eq. 
98]), we get 



R(k,u) = 



(2ir) d u{k)k 2 p 



{iku Q + uj) 2 - {2i:) d u(k)k 2 p 



(103) 



The response of the system to a pulse is given by Eq. (|14p with Eq. (|103|) . Taking its inverse Laplace transform, we 
obtain 



6<f>{k,t) 



(2ir) d u(k)k 2 p 



lc 27r (w + iu Q k + v / (2TT) d u(k)k 2 p)(uj + iu k ~ y/ ' (2n) d u(k)k 2 p) ' 
This integral can be easily performed with the residue theorem. For a repulsive potential (u(k) > 0), 

6®(k,t) = -^(2ir) d u(k)k 2 psm (j ' (2ir) d u(k)k 2 p?J e~ Uokt , 

and for an attractive potential (u(k) < 0), 



8$(k,t) = ^(2Tr) d \u(k)\k 2 psmh \^J{27r) d \u(k)\k 2 ptj e - u ° kt . 
The evolution of the perturbation is consistent with the discussion given in Section IVAl 



(104) 



(105) 



(106) 



D. The response to a step function 



The response of the system to a step function is given by Eq. (fl9|) with Eq. (|103[) . Taking its inverse Laplace 
transform, we obtain 



(2iT) d u(k)k 2 p 



/c 2tt u (uj + iu k + v / {2Tr) d u(k)k 2 p){oj + iu Q k - ^ (2ir) d u(k)k 2 p) 
This integral can be easily performed with the residue theorem. For a repulsive potential {u(k) > 0), 

^(k,i)=<5$oo(k)Jl-e- llofct 



(2Tt) d u(k)p 
and for an attractive potential (u(k) < 0), 



sin J(2Tr) d u(k)k 2 pt + cos J (2it) d u(k)k 2 pt 



5*(k,t) = 5*oo(k) |l- e -" ofc * 
where 



sinh J{2ir) d \u{k)\k 2 pt) + cosh J (2ir) d \u(k)\k 2 pt 



(2ir) d \u(k)\p' 



S3, (kl (2,) d u{k)pk 2 . 



(107) 



(108) 



(109) 



(110) 



In the stable case, the zeros of the dielectric function lie in the lower half-plane and the perturbation asymptotically 
tends towards the distribution defined by Eq. (|110[) . 
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VI. LINEARIZED HYDRODYNAMIC EQUATIONS 

In this section, we derive hydrodynamic equations satisfied by the perturbed density Sp(k,t) in the linear regime. 
We consider the case of a collision-dominated gas (Euler) and the case of a collisionless system (Vlasov) and show the 
correspondence between these two systems. 



A. Linearized hydrodynamic equations for a collision-dominated gas 

For a collision-dominated gas, it is straightforward to obtain the equation for the perturbed density 6p(r, t). Taking 
the time derivative of Eq. (|29[) and the divergence of Eq. p0[) . and combining the resulting equations, we obtain 



d 2 Sp dSp 



dt 2 



f - cjASp = P AS<Z> + pA$ e 



(111) 



This equation without the right hand side is called the telegraph equation. However, Eq. (|111[) is more complicated 
than the telegraph equation because 5$ depends on Sp according to Eq. (|3"Tj) . It is therefore an integro-diffcrcntial 
equation. Taking the Fourier transform of Eq. (|lllj) . we get 



where <5$ is related to Sp by Eq. (|34[). The foregoing equation can be rewritten as 

d 2 Sp t dSp 2 



dt 2 ' ■f^ 1 + w o(*0^ = -pfc a *e. 



For £ = (Euler), it reduces to 



For £ — ^ +oo (Smoluchowski), we get 



d 2 5p 



+ ujl(k)Sp = -pk 2 $ e . 



^+^ 2 Q (k)Sp=-pk 2 ^ 



(112) 



(113) 



(114) 



(115) 



B. Linearized hydrodynamic equations for a collisionless system 

For a collisionless system, the perturbed density Sp(k, t) satisfies the integral equation [see Eq. (|89[) ]: 

Sp{k,t)=iJ dt / [<y$(k,0 + *e(k I t / )] Jdvk- l^e"*^-*'), 
where 5$(k, t') is related to Sp(k,t') by Eq. ([74")) . Integrating the last term by parts, we obtain 



Sp{k, t) = -k 2 J dt' |£$(k, t') + $ e (k, i')J (t - / rfv /(v)e 



-ik-v(t-t') 



Introducing the reduced distribution function f(v), we get 



Sp(k,t) = -k 2 dt' \s$(k,t') + $e(M')l (i-i') / dw/(w)e- <fco (*-*'>. 



e- 4 «"/(t;), 



Defining the Fourier transform of f(v) in velocity space by 



/(<?) 



2tt 



(116) 



(117) 



(118) 



(119) 
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we can write Eq. (|1 18[) in the form 



Sp(k,t) = -2nk / dt' <5$(k,i') + $ e (M') qf{q), 



where q = kit — t'). Taking the first and second derivatives of this equation with respect to time, we get 



d5p 



= -2-nV 



dt' 



5$(k,t') + $ e (k,t')\ ^kf(q)l 



(120) 



(121) 



d 2 6p 



-pk 2 



<5$(k,i) + $ e (k,t) 



dt' 



5*(k,t / ) + *e(k,t / ) 



dq- 



■[9/(9)1 



(122) 



To obtain Eq. (|122j) . we have used the fact that (d/dq)[qf(q)} = qf'(q) + f(q) is equal to /(0) = p/2n when g = 0. 
In general, we cannot obtain a closed partial differential equation for Sp(k, t), unlike in a collision-dominated gas. An 
exception concerns the waterbag and Cauchy distributions. 



1. The waterbag distribution 

Taking the Fourier transform of the waterbag distribution defined in Section [IV El we obtain 

/(9) = f^^. (123) 
27r v m q 



Substituting the identity 
in Eq. ([122)) . and using Eq. (|T20| . we get 

d 2 Sp 

Comparing this equation with Eq. (|112[) . we see that a collisionlcss system with the waterbag distribution behaves 
as a collision-dominated gas in which the velocity of sound c s is replaced by the maximum velocity v m . This system 
does not experience Landau damping. We note that the dispersion relation (|94[) of the waterbag distribution can be 
obtained directly from the hydrodynamic equation (|125p . 



— [ q f(q)]=-vl q f( q ). (124) 



v 2 m k 2 6p= -pfc 2 (5$ + $ e ). (125) 



2. The Cauchy distribution 
Taking the Fourier transform of the Cauchy distribution (|96p . we obtain 

2tt 



f(q) = J?-e- u °«. (126) 



2 [qf(q)} = -2u -[qf(q)} - u 2 qf(q), (127) 



Substituting the identity 

jL [q f( q )] = -2u f q 
in Eq. |[T2^|) . and using Eqs. p^O)) and CCD]) , we obtain 

^ + 2u k^- + u 2 k 2 Sp = -pk 2 (5$ + $ e ). (128) 

Comparing this equation with Eq. (|112p . we see that a collisionless system with the Cauchy distribution behaves as 
a collision-dominated gas in which the velocity of sound c s is replaced by u$. Furthermore, this equation exhibits a 
damping in Fourier space with a damping coefficient proportional to k. This is of course a manifestation of Landau 
damping. However, it appears here in a simple manner without having to carry out contour integration in the 
complex plane. We note that the dispersion relation (|99[) of the Cauchy distribution can be obtained directly from 
the hydrodynamic equation (|128[) . 
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VII. ASYMPTOTIC VALUE OF THE PERTURBATION 



The response of the system to a step function is given by the general expression (|21[) . In the case where all the 
zeros of the dielectric function e(k, to) strictly lie in the lower half plane (a;, < 0), an assumption that will be made 
in this section, the perturbation <5$(k, t) relaxes towards the asymptotic value <5$oo(k) given by Eq. (|22|) . It is clear 
that this asymptotic value corresponds to the steady state of the system under a weak external potential $ e (k) 1 . We 
can make this expression more explicit for a collision-dominated gas and a collisionless system. 



A. Collision-dominated gas 

For a collision-dominated gas, the dielectric function (j38|) with uj = reduces to 

e(fc> 0) = l + M^£ . (129) 

Substituting this expression in Eq. (|22|) . we get 

c- „s (2n) d u(k)pk 2 ~ s , . 

We can show that Eq. fjl 30[) corresponds to the steady state of the damped Euler equations ([23]) - ([251) with the 
equation of state p(p) under a weak external field $ e (r). Indeed, in the presence of an external field, the condition of 
hydrostatic equilibrium reads 

Vp + /?V$ + pV$ e = 0. (131) 

If the applied field is weak, we can linearize Eq. (|131[) about the unperturbed distribution. If the unperturbed system 
is spatially homogeneous, we obtain 

cJVtfp + pVSS + pV* e = 0. (132) 

After integration, we get 

c 2 Sp(r) + p(M>(r) - (5$)) + p($ e (r) - (* e » = 0, (133) 

where the constant of integration has been determined by using the fact that the external potential does not change the 
mass (the brackets denote a space average over the whole domain) . Taking the Fourier transform of this expression, 
and using Eq. (|34|) . we get Eq. (|130[) . Therefore, when submitted to a step function, a stable collision-dominated 
gas converges for t — > +oo towards an asymptotic distribution (5$oo(k) which is the steady state of the damped Euler 
equations (|23l) -([25 |) under a weak external field <E> e (r) 2 . 

B. Collisionless system 

For a collisionless system, the dielectric function ([75} with ui = reduces to 

e(fc, 0) = 1 - (27t) d u(k) J ffi^ dv. (134) 
Substituting this expression in Eq. (|2"2"]l , we get 

1 - {2n) d u{k) J J-jp- dv 



1 It suffices to take d/dt = in the linearized equations of Sections IIII Al and IIV Al The equations of this section remain valid provided 
that we use just Fourier transforms and set ui = 0. Eqs. J35 D and H75H then yield Eq. (1221) . 

2 This is true for the damped Euler equation (§ > 0) since all the zeros of the dielectric function lie in the lower half-plane. For the Euler 
equation (£ = 0), the perturbation oscillates indefinitely about ^^(k) as explained in Section IlII D II 
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When /(v) = F(v 2 /2), we can show that Eq. (|135[) corresponds to the steady state of the Vlasov equation (|69[) - 
([70)) with the distribution F under a weak external field 3> e (r). Indeed, from the spatially homogeneous distribution 
function /(v) = F(v 2 /2), we can define a spatially inhomogeneous distribution function /(r, v) = F(e) with e = 
v 2 /2 + $*(r) where $t(r) = 3?(r) + ^(r) is the total potential. This distribution /(r, v) = F(e) is a steady state 
of the Vlasov equation. We have to show that, when the external potential is weak, this distribution function leads 
to Eq. (|135[) . To that purpose, we first note that, for any collisionless system with a distribution function of the 
form / = F(e), there exist a corresponding barotropic gas p = p(p) with the same density profile (see Appendix [A"[) . 
Furthermore, the condition / = F(e) implies the condition of hydrostatic equilibrium (|131[) . Therefore, when the 
external potential is weak, proceeding as in Section [VII Al we find that the perturbation is given by 



Now, the velocity of sound in the corresponding barotropic gas is given by Eq. (|A3[) . Substituting this relation in Eq. 
(11361) . we recover Eq. (|135p 3 . Therefore, when submitted to a step function, a stable collisionless system converges 
for t — > +oo towards an asymptotic distribution (5$ OC) (fc) which is the steady state of the Vlasov equation (|69|) - (f7U|) 
with the distribution F under a weak external field <I> e (r) 4 . Using Eq. (| A3|) . we note that the dielectric function f| 134[) 



can be written as 



just like in a collision-dominated gas. 



e(M)=1+ e«, m 



C. The relation with the correlation function 



For the isothermal distribution function (|145[) . corresponding to statistical equilibrium, the asymptotic expression 
of the perturbation reduces to 

*-« " - t£%X)», *-U- < 138 > 

As a confirmation of the previous calculations, we can directly establish that (5$oo(k) corresponds to the potential 
produced by the Boltzmann distribution with a weak external field. For a system at statistical equilibrium under an 
external field, the Boltzmann distribution function reads 

/(r,v) = ^ e -^[- 2 /2+*W+*a(r)]. (139) 

Integrating over the velocity, we get 

e -/3[*(r)+* e (r)] 

p{ r )= N ~r- rnr-rx\ ; ■ (!40) 

If the external field is weak, we can linearize the previous equation. Assuming that the unperturbed system is spatially 
homogeneous, we obtain 

Sp(r) = -0p(*$(r) - (5$)) - /3p($ e (r) - (* e ». (141) 

This returns Eq. (|132[) with c 2 = T, leading to Eq. (|138[) . Therefore, Eq. f|138[) corresponds to the statistical 
equilibrium state of the system under a weak external potential. We stress, however, that Eq. (|135[) [rcsp. Eq. (|136[) ] 
is more general since it is valid for an arbitrary distribution f(v) [resp. /(v) = F(v 2 /2)], not only for the Boltzmann 
distribution. 



3 If we expand F(e) = F(v 2 /2 + &t{r)) for &t(r) 1, we obtain Eq. H76H with oj = 0. After integration over v, this leads to Eq. 11751 
which is equivalent to Eq. II135I I. This shows the compatibility of the different approaches. 

4 This is true for most distribution functions (that experience Landau damping). A notorious exception is the water bag distribution that 
does not experience Landau damping [T3| . In that case, the perturbation oscillates indefinitely about 5&oa (k) as in a perfect gas (£ = 0). 
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Finally, if we recall the expression of the Fourier transform of the two-body correlation function of a spatially 
homogeneous system with long-range interactions [24], Ifi3| : 

<0^r,h(h\ (27r)V(fc)fc 2 l-e(fc,0) 

(2 " } nh(k) - Tk 2 + (2n) d u(k)k 2 p = e(k,0) ' (M2) 
and compare its expression with Eq. (|138|) . we obtain the relation 

<$$oo(k) = (2ir) d h{k)n$ e {k), (143) 

where n = p/m in the numerical density. Therefore, the measure of the asymptotic potential produced by a system 
submitted to a step function allows us to determine the equilibrium correlation function of the system. 



VIII. APPLICATION TO SPECIFIC DISTRIBUTION FUNCTIONS 



In the previous sections, we have given general criteria of stability and general expressions for the asymptotic distri- 
bution of a stable homogeneous system submitted to step function. These results are valid for arbitrary distributions 
and for arbitrary potentials of interaction. Furthermore, when /(v) = F(v 2 /2), we have shown that the results take 
the same form in a collisionless system and in the corresponding barotropic gas. In particular, they can be expressed 
very simply in terms of the velocity of sound [see Eqs. (j53"j) and (|136[) ]. As explained in Appendix \K\ the velocity of 
sound is a function of a (generalized) temperature T. In turn, the temperature can be expressed in terms of the energy 
E by Eq. (|A4|) . In this section, we illustrate these results for specific distribution functions (isothermal, polytropic, 
watcrbag, and Fermi-Dirac). In the context of the HMF model, we can derive simple formulae that complete those 
obtained previously in Refs. flH [ill]. In particular, we obtain the general result 

M °° = 2^TT' c * > 5> (144) 

for the asymptotic magnetization 5 . The asymptotic magnetization M m is proportional to the magnetic field h, and 
the magnetic susceptibility x = M^jh is given by a generalized Curie- Weiss law: \ = l/(2<^ — 1). 



A. Isothermal distribution 



For the isothermal distribution function 

d/2 



/( v ) = (0' pz-^ v2 > 



we obtain the equation of state 

P = pT, (146) 

and the velocity of sound 

c 2 s = T. (147) 

For the HMF model, we can substitute this result in Eq. (| 144|) to express the stability criterion and the asymptotic 
magnetization in terms of the temperature T. We obtain 

M --W^f:y T>Tc 4 (148) 



5 For the HMF model, it is more convenient to work in terms of the magnetization M than in terms of the potential <5<£>(0). They are 
related to each other by 5^(9) = —McosB yielding <5<i> n = — (M/2)8 n ,±l- It is also convenient to write the external potential in the 
form <3? e = —hcosB, yielding (4 e ) n = —(h/2)S n ,±i, where h can be interpreted as a "magnetic" field. 
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where we have used E = T/2 + 1/2 [see Eq. ([Al])]. This returns the results obtained in [l^, [l^] . Eqs. (|TI5 )) - p^|) 
are also identical to the expressions of the magnetization of the HMF model at statistical equilibrium under a weak 
magnetic field h (see Eqs. (100) and (107) of |6J])- As we have explained in Section fVIIi this property is general. 

B. Polytropic distributions 

The polytropic distributions can be written as [65j : 



/4 



1/n _ v 2 /2 



{n + l)K 



n-d/2 

(150) 



where Z is given for n > d/2 by 



z = s ^-. rW2)r(i-y + n) 

1 (n + 1) 



and for n < —1 by 



Z = S d 2 d /^ T ^ {d,2 h -K(n + l)] d '\ (152) 
1 (a/2 — n) 

The constant K is called the polytropic temperature. These distribution functions were also introduced by Tsallis 
[6r| in his generalized thermodynamics. The corresponding equation of state is 

y = Kp 1 1 7 = 1 + -. (153) 
n 

The velocity of sound is given by 

c 2 s = Kryp"*- 1 . (154) 

For the HMF model, we can substitute this result in Eq. (| 144|) to express the stability criterion and the asymptotic 
magnetization in terms of the polytropic temperature K. Since p = 1/2tt in the homogeneous phase, it is convenient 
to define the polytropic temperature by & = K{l/2ir) 1 ~ 1 so that c 2 s = 70. Then, we obtain 

M —2^~V (155) 

= 47^27-1 ' E> ^ +1 r ^ 

where we have used E = 0/2+1/2 [see Eq. (|A4[) ]. For n — s- +oo, the distribution function (|150|) reduces to the 
isothermal distribution (|145[) and we recover the results of Section IVIII Al 

C. Waterbag distribution 

The waterbag distribution, defined by / = 770 for — v m < v < v m and / = otherwise, is a particular polytrope of 
index n — d/2. The corresponding density and pressure arc p = rjoSdvf n /d and p = rjQSdvf^ 2 /[d(d + 2)], leading to 
the polytropic equation of state 

d + 2 1 / d \ 2/d , , 

,-K,-, 7 = — K = — ( — ) . (157) 
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The velocity of sound is c 2 s = vf n /d. In d = 1, we have c s = v m in agreement with the results of Section llVEl For the 
HMF model, we can substitute these results in the general equation (|144[) to express the stability criterion and the 
asymptotic magnetization in terms of v m . We obtain 

M - = v - > 71' (158) 



= E > ^ (159) 

,/6 + 1/2 [see Eq. (|A4[) ]. We recall, however, that a collisionless system with the waterbag 
distribution submitted to a step function does not relax towards the magnetization M m but oscillates about it as 
M(t) = Moo[l - cos(w i)]with the pulsation uj = ±^Jv 2 ~\J 7 l (see Eq. (J54J) with the remark of Section HVE| . This 
result was obtained in [l5| and it was confirmed by direct numerical simulations. When submitted to a pulse, the 
magnetization behaves as M(t) = ~ sin(u)ot)/u>o according to Eq. f|44|) . 



D. Fermi-Dirac and Lynden-Bell distributions 

For the Fermi-Dirac distribution function (see, e.g., 



/w - ttS^' (160) 



we obtain 



P= ^/2-l( A ). P= df3d/2+1 Id/»W, (161) 

where rjo is the maximum value of the distribution function fixed by the Pauli exclusion principle, S d is the surface of 
a unit sphere in d dimensions, and I n (X) = J + °° + Xe x ) dx are the Fermi integrals. This distribution function 

also corresponds to the prediction of Lynden-Bell [68[ in his statistical theory of the violent relaxation of the Vlasov 
equation. In that context, rjo represents the initial value of the distribution function (in the two-levels case). Equations 
(11601) and (|161|) define the equation of state p(p) of the Fermi gas in parametric form. In the non degenerate limit 
(T — > +00), the Fermi-Dirac distribution function reduces to the classical isothermal distribution function (|145|) . The 
corresponding equation of state is p — pT (see Section TVIII A|) . In the completely degenerate limit (T —> 0), the 
Fermi-Dirac distribution function reduces to the waterbag distribution / = rjo for — vp < v < vf and / = otherwise, 
where vf is the Fermi velocity (see Section IVlII C|) . The corresponding equation of state is the polytropic equation of 
state l|157p with v m = Vp- The velocity of sound can be written as 

For a given temperature T, we can determine A from Eq. (|16H -a) and c s from Eq. (|162[) . Therefore, the velocity 
of sound is a function of the temperature: c s = c s (T). In the non degenerate limit, c 2 s = T and in the completely 
degenerate limit c 2 s = vp/d. For the HMF model, we can substitute this result in Eq. (|144[) to express the stability 
criterion and the asymptotic magnetization in terms of the temperature T. We obtain 

M ~ = 2^FT< ^ T)> \- (163) 

Using Eqs. (|A4j) and (|161I b). we find that the energy is given by 

Eliminating A between Eqs. (|164p and (|161l -a). we can obtain E(T) in parametric form. We can then express the results 
of Eq. (|163|) in terms of E instead of T. In the Lynden-Bell theory, 770 is an additional control parameter related to the 
initial condition, so that the results actually depend on (T, 770) or on [E, t)q). The stability of the spatially homogeneous 
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distribution [see criterion (|163l -b)] was studied in [30j | . For T = 0, corresponding to E = E m i n = 1/(967t 2 ?7q) + 1/2, 
we find that the system is stable for 770 < l/(2v27r). In that case, = h/[l/ (Sir 2 ^) — 1] (this corresponds to the 
results of Section fVIII CI with v m — vp = l/(47r?yo)). The asymptotic expression of the magnetization M^iT) in the 
limit T — > or E — > E min (for fixed rj ) was obtained in [l5j . Using the Sommerfeld expansions of the Fermi integrals 
(see, e.g., Eqs. (39) and (40) of |13|), we can check that Eq. (|163[) returns the results of [l5|. However, Eq. f)163|) is 
more general, as it is valid at any temperature T > or at any energy E > E m i n . 



IX. THE INITIAL VALUE PROBLEM 



In Section IIV1 we have studied the response of a collisionless system described by the Vlasov equation to a weak 
external potential $ e (r, i) using the linear response theory. Here, we compare these results with those obtained when 
the system is isolated (i.e., $ e (r, t) = 0), but the distribution function is slightly perturbed at t = 0. This is the 
so-called initial value problem of the linearized Vlasov equation Q ■ 

Taking the Fourier-Laplace transform of the linearized Vlasov equation (|7l"j) - (|72p . and assuming now that 
Sf(r, v, 0) 0, we obtain 

*/(k,v >W ) = J^S-fifr*;) + .ffi' V,0 \ , (165) 
k • v uj i(k ■ v — uj) 

where Sf(k, v, 0) is the Fourier transform of the initial perturbation Sf(r, v, 0). Integrating Eq. (|165[) over the velocity 
and using Eq. ([74)) . we find that 

^(k,.).(2.)^/|^dv. (166) 
e(k, uj) J i{k-v — w) 

Substituting this result back into Eq. (|165[) . we get 

v ' k-v-w v ; e(k,w) J i(k-v'-cj) i(k • v - u) 

This is the exact solution of the initial value problem for the linearized Vlasov equation. It is conveniently written in 
terms of a resolvent operator that connects 5/(k, v, ui) to the initial value 

Sf{k,v,u)) = y'(iv'i?(v|v',k,w) ( 5/(k,v,0), (168) 

with 

W,k,.) = JLg- (2 ^ «( fe ) \ + . (169) 

k • v — to e(k, uj) i[k -v — uj) i(k ■ v — uj) 
If we consider an initial condition of the form Sf(r, v, 0) = m6(v — v')<5(r — r'), leading to 

5f(k,v,0) = j^S(v- v')e- iW , (170) 



wc find that 



®j_ L,\ n-ik-r' ™ „— ik r' 



g/(k,v,a,)= k '«v y m 6 7' +^( V - V /) « (171) 

k-v — we(k, w) z(k ■ v' — uj) (27r) ct z(k ■ v — uj) 

On the other hand, if we substitute Eq. (f74|) in Eq. (|165|) . wc obtain 

Sf(k, v, «) = , k '^ (2n) d u(k) { Sf(k, v', uj) dV + .^ (k ' V ' 0) . . (172) 
k-v — uj J i[k ■ v — uj) 

This is an integral equation for <5/(k, v,w) that is equivalent to Eq. (|167[) . 

These results are well-known in plasma physics ji|. They have been recalled in order to facilitate the comparison 
with the results obtained with the linear response theory. In particular, Eq. (|166p may be compared to Eq. (|75[) and 
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Eq. (|167j) may be compared to Eq. (|76l) . In general, the two sets of equations differ, showing that the initial value 
problem is not equivalent to the linear response theory. However, if we consider an initial disturbance of the form 

<y/(k,v,0)=ik-|^$ e (k) ) (173) 

where $ e (k) is an arbitrary function of k (independent on v), we see that Eqs. (|73|) and (|165|) coincide provided that 
$ e (k, w) = $ e (k). Since the Laplace transform of the external potential is independent on u, it corresponds to a 
pulse with an amplitude $ e (k) (see Section Til Bp . Therefore, the effect of an initial disturbance of the form (|173[) on 
an isolated system is equivalent to submitting this system to a pulse $ e (k, t) = $ e (k)<5(£). 

The perturbation <5$(k, uS) given by Eq. (|166|) appears as a product of two factors: e(k, and an integral 
involving the initial condition <5/(k, v, 0). The first factor is due to collective effects, and its role in the evolution 
of the perturbation has been already discussed in Section III Bl This term can produce damped, steady, or growing 
oscillations. On the other hand, the integral corresponds to the excess density produced by an initial disturbance in 
a gas of non-interacting particles (i.e., for which e(k, oS) = 1). It is typical of an individual particle behavior. We 
expect that the effect of this term will disappear for late times. Actually, it can be shown that the integral produces 
damped oscillations (i.e. its poles are necessarily in the lower half plane or on the real axis) This is associated 
with the phenomenon of "phase mixing" which is an irreversible homogenization process in which the interactions 
play no role 6 . In a sense, this term is the counterpart of the polarization function in the linear response theory. 

The individual particle behavior corresponds to the integral in Eq. (|166[) . This integral represents the excess density 
in the absence of interaction. It can be written as 

f 6f(k,v,0) J 

5p {k,w)= / — rcif, (174) 

J L i(kv-uj) 

where L denotes the Landau contour and Sf(k,v,0) is the reduced distribution function corresponding to the initial 
disturbance (i.e. it has been integrated over the coordinates perpendicular to k • v) . The individual particle behavior 
can be studied analytically for an initial disturbance of the form 

A(k) 1 

6f(k,v,0) = — ^ , rg-. (175) 

w l 

This example has been discussed in Ref. Using the residue theorem, we easily obtain 

5po(k,uj) = — A(fc) (176) 

1(KVq — IKWq — 10) 

Then, taking the inverse Laplace transform of this expression, we get 

Sp (k, t) = A(k)e- ikvot e- kwot . (177) 

This corresponds to a damped oscillation: The perturbation propagates with a group velocity vq and dies out ex- 
ponentially rapidly on a timescale = (fcwo) -1 depending on the wavenumber k. We now treat another example 
that has not been discussed before (to the best of our knowledge). We consider an initial disturbance of the form 
6f(k,v,0) = A(k)/2w m if — w m < v < w m and 5f(k,v,Q) = otherwise. The excess density (|174[) can be written 
explicitly as 



2w m ik \u> + kw 



If we consider the perturbed distribution function I I167I I. we sec that there is a real pole u) = k ■ v. It produces an undamped oscillation 
exp(— ik ■ vt) whose pulsation is proportional to the velocity v of the particles. It is therefore an individual particle effect. Let 
us consider the excess density <5/5(k, t) = J <5/(k, v, t) dv. The contribution to this integral from the various velocities will produce 
destructive interferences of the oscillations, hence this part of <5p(k,t) tends to zero for long times: this is the phenomenon of phase 
mixing. The other poles of Eq. I|167|l correspond to the zeros of the dielectric function: e(k, cj) = 0. They depend on k but not on v. 
They describe the collective behavior of the system. They produced damped or growing oscillations that "resist" the integration over v 
1- 
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Taking the inverse Laplace transform of this expression, and integrating by parts, we obtain 

sm(kw m t) 



Sp (k,t) = A(k)- 



kw m t 



(179) 



We see that the damping of the initial perturbation is algebraic, behaving as t^ 1 for large times, instead of being 
exponential as in Eq. (|177p . It also exhibits oscillations with a pulsation kw m that depends on the wavenumber k. 
As a result of this slow damping, the initial disturbance will have a long-term effect. 

It is instructive to recover these results in a different manner Q . Taking the Laplace transform of Eq. f|l T4[) before 
performing the velocity integration, we get 



Sp (k,t) = I 5f{k,v,0)e- ihvt dv 



(180) 



Under this form, we see that the evolution corresponding to the free-particle motion consists of a superposition of 
waves exp(— ikvt) propagating with a group velocity equal to the velocity of the individual particles. These waves 
will in general interfere destructively as time goes on, leading to the free-motion damping discussed above. Actually, 
the integral (|180[) can be calculated easily returning the results (|177|) and (|179[) . 

Finally, as shown in Q , when the unperturbed distribution function is the Cauchy distribution (|96[) , and when the 
initial disturbance in given by Eq. (|175|) . it is possible to determine the evolution of the perturbation 5&(k,t) given 
by Eq. (|166p analytically. Generalizing the calculation to an arbitrary potential, we get 



6$(k,t) = (2w) d u{k)A(k)\ e -^v a t-kw Q t 



Q, 2 



(kvo + i(uo — wo)k) 2 — O 2 
kvQ + ik(uo — wq) — Q kvo + ik(ua — wq) + Q 



an 



(181) 



where we have defined the pulsation £l(k) = \J (2ir) d u(k)k 2 p which is real for a repulsive interaction (u(k) > 0) and 
purely imaginary for an attractive interaction (u(fc) < 0). The first term in braces in Eq. (|181[) corresponds to 
individual effects and the second term corresponds to collective effects. 



X. CONCLUSION 



We have applied the linear response theory to systems with long-range interactions. Our study completes previous 
investigations [HI, [l6|. We have considered a collision-dominated gas (Euler) and a collisionless system (Vlasov). We 
have shown that the response of these systems to an external field is in general different except for the waterbag 
distribution. In that case, there is no Landau damping and the dielectric function of the waterbag distribution 
coincides with the dielectric function of a collision-dominated gas without dissipation (£ = 0). When submitted to a 
step function, these systems oscillate permanently without reaching a steady state. For more generic distributions, 
such as the Cauchy distribution, there is Landau damping and the evolution of a collisionless system described by 
the Vlasov equation resembles the evolution of a collision-dominated gas described by the damped Euler equation 
(£ > 0). When submitted to a step function, these systems relax towards a steady distribution. When / = F(v 2 /2) 
this distribution is the same in the collisionless system and in the corresponding barotropic gas. It corresponds to the 
steady state of the system under a weak external field. However, the relaxation towards this steady state is in general 
different in a collisional gas and in a collisionless gas because their dielectric functions differ. We have also considered 
unstable systems. In that case, the linear response theory is valid only for short times, before the perturbation has 
significatively grown. Physical applications of the linear response theory will be given in a future paper [69^ . 



Appendix A: The corresponding barotropic gas 

We consider a collisionless system described by the Vlasov equation (|69 |) -([70 |) . We assume that the system is 
spatially homogeneous. Any distribution function /(v) is a steady state of the Vlasov equation. We restrict ourselves 
to distribution functions of the form /(v) = F(v 2 /2). From the spatially homogeneous distribution function /(v) = 
F(v 2 /2), we can define a spatially inhomogeneous distribution function /(r,v) = F(e) where e = v 2 /2 + <I>(r) is the 
individual energy of the particles. This distribution /(r, v) = F(e) is also a steady state of the Vlasov equation. 

For any collisionless system with a distribution function of the form /(r,v) = F(e), there exist a corresponding 
barotropic gas with the same equilibrium density (70j . Indeed, introducing the density p(r) = f F(e)dv and the 
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pressure p(r) = \ J F(e)v 2 dv, we have p(r) = p[<f>(r)] and p(r) = p[<f>(r)]. Eliminating the potential $(r) between 
these two expressions, we obtain a barotropic equation of state p = p{p) that is entirely determined by the function 
F. We can easily show that the condition /(r, v) = F(e) implies the condition of hydrostatic equilibrium. Indeed 



Vp = ~ V$ J F'{e)v 2 dv = iv$ J |^ • v dv = -V$ / / dv = -pV$. 



(Al) 



Finally, we can relate the velocity of sound c 2 = p'(p) in the corresponding barotropic gas to the distribution function 
/(v). To that purpose, we first note that the condition of hydrostatic equilibrium (|A1[) combined with the equation 
of state p = p{p) implies that = —/?($). Therefore, 
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^ /^(e)dv /ff(«,*)i(fo (A2) 

In the last expression, we have integrated over v\...Vd-i and noted v for vg, and f(v,&) for J f dvi...dvd-i- For a 
spatially homogeneous system, the velocity of sound takes the simple form |2j |22| : 

c * = "r+»m du ' (A3) 

J — oo V 

where f(v) is the reduced distribution of Section HV Al 

In general, the distribution function f(v, p, T) and the pressure p(p, T) depend on the density p and on an external 
parameter T which can be identified with a (generalized) temperature. Therefore, the velocity of sound is a function 
c s (p,T) of the density and of the temperature. Substituting this relation in the general equations ([83)) and (|136[) . we 
can express the stability criterion and the asymptotic distribution in terms of T and p. On the other hand, the total 
energy of a spatially homogeneous system can be written as 

E =\J f v ' 2 drdv + W = ^ J P dr + W = ^ P (p, T)V + Wo, (A4) 

where V is the volume and Wo the potential energy in the homogeneous phase. From this equation, we can relate the 
energy to the temperature and express the results (|83p and (|136p in terms of E. Some explicit examples are given in 
Section EUU 



Appendix B: Plasmas, self-gravitating systems, and the HMF model 



In this Appendix, we provide explicit solutions of the dispersion relation for plasmas, self-gravitating systems, and 
for the HMF model. We consider the waterbag and the Cauchy distributions. 



1. Coulombian plasmas 



For a 3D Coulombian plasma, the potential of interaction is the solution of the Poisson equation Ait = 
— 47r(e 2 /m 2 )<5(x) yielding (27r) 3 w(fc) = 4ne 2 /m 2 k 2 . We introduce the plasma pulsation wp = (Aitpe 2 /m 2 ) 1 / 2 . 

For the Cauchy distribution, the dispersion relation becomes lu = ±u>p — iuok. The perturbation oscillates 
with a pulsation uj r = ±tup and is damped at a rate uii = —uok < 0. For the long wavelengths, the damping is 
negligible and the plasma oscillates with the proper pulsation tup. This reflects the collective behavior of the system. 
The natural limit of collective behavior corresponds to up ~ u k leading to the Debye wavenumber kjj = lop/uq = 
(4Trpe 2 /m 2 u 2 ) 1 / 2 . 

For the waterbag distribution, the dispersion relation (|94|) becomes to 2 = v^k 2 +wp. The perturbation oscillates 
with a pulsation ui = ±yjv 2 n k 2 + ulp without attenuation. For the long wavelengths, the plasma oscillates with the 
proper pulsation u)p. The Debye wavenumber is kp, = iop/v m = (4:7rpe 2 /m 2 v 2 n ) 1 ^ 2 . 



2. Self-gravitating systems 



For a 3D self- gravitating system, the potential of interaction is the solution of the Poisson equation Au = 4irGS(x) 
yielding (2*7r) 3 u(fc) = —AirG/k 2 . We introduce the gravitational "pulsation" luq = ^AnGp. 
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For the Cauchy distribution, the dispersion relation (|99|) becomes u> = ±zwg — iuok. The Jeans wavenumber is 
kj = uj g /u = (AirGp/ul) 1 / 2 . The syste m is stable for fc > fcj and unstable otherwise. In the stable case, the 
perturbation is damped with an exponential rate u>i — ±wg — u^k < 0. In the unstable case, the perturbation grows 
with an exponential rate uii = uiG — Utyk > (the other mode is damped with an exponential rate uii = —uiQ — uok < 0). 

For the waterbag distribution, the dispersion relation ([94]) becomes u> 2 = v 2 n k 2 — u>q. The Jeans wavenumber is 
kj = Lua/vm = (AirGp/v 2 ^ 1 / 2 . The system is stable for k > kj and unstable otherwise. In the stable case, the 
perturbation oscillates with a pulsation oj r = ±^Jv 2 n k 2 — ui 2 G . In the unstable case, the perturbation grows with an 
exponential rate uii = y'ug — v 2 n k 2 > (the other mode is damped with an exponential rate Wj = — \/Uq — v^k 2 < 0) . 



3. The repulsive HMF model 

For the repulsive HMF model, the potential of interaction is u = jp [1 + cos(# — 8')] yielding u n = ^(2S n! o + <J n ,i + 
<5„._i). In a collisionless system (Vlasov), only the modes n — ±1 can propagate [22l |. 

For the Cauchy distribution, the dispersion relation (|99|) becomes ui = ±l/v2 — iuq. The perturbation oscillates 
with a pulsation ui r = ±l/v7 and is damped at a rate u>i = — uo < 0. For uq — > 0, the damping is negligible and the 
system oscillates with the proper pulsation l/\/2. This reflects the collective behavior of the system. The natural 
limit of collective behavior corresponds to Uq = 1/2. 

For the waterbag distribution, the dispersion relation ([94]) becomes to 2 = + 1/2. The perturbation oscillates 
with a pulsation uj = ± v 2 n + 1/2 without attenuation. For v m —} 0, the system oscillates with the proper pulsation 
1/V2. The natural limit of collective behavior corresponds to v 2 n = 1/2. 



4. The attractive HMF model 

For the attractive HMF model, the potential of interaction is u = i [1 — cos^ — 9')] yielding u n = j^(25 n fi — 
S n ,i — S n ,-i)- In a collisionless system (Vlasov), only the modes n = ±1 can propagate |22l |. 

For the Cauchy distribution, the dispersion relation (|99|) becomes w = ±i/ V2—iu n . The system is stable for > 1/2 
and unstable otherwise. In the stable case, the perturbation is damped with an exponential rate uii = — < 0. 

In the unstable case, the perturbation grows with an exponential rate uii = l/s/2 — «o > (the other mode is damped 
with an exponential rate u>i = — 1/\/2 — ito < 0). 

For the waterbag distribution, the dispersion relation (|94[) becomes uj 2 = v 2 n — 1/2. The system is stable for 
v m > l/v2 and unstable otherwise. In the stable case, the perturbation oscillates with a pulsation uj r = i-y/vj^ — 1/2. 
In the unstable case, the perturbation grows with an exponential rate uii = -\/l/2 — > (the other mode is damped 
with an exponential rate u>i = — y/ 1 / 2 — v 2 n < 0). 



Appendix C: An identity 

For a collisionless system described by the Vlasov equation, the Fourier-Laplace transform of the correlation function 
of the fluctuations of the potential is given by Eq. (28) of Ref. [tTJ] . Taking the inverse Laplace transform of this 
equation and using Eq. (j74)) . we obtain 

(^(MMk',t')> = ^(k + k') / «- flt - v(t - | e(k ^ v) v) , a dv- (ci) 

The equilibrium correlation function is therefore 

{SffkWW) - + k') / , e(k ^ V) v) | 2 dv, (C2) 

where /(v) is the Maxwell distribution (| 145[) . On the other hand, the correlation function can be written as (see, e.g., 
Appendix A of [Uj]): 

(8p(k)Sp(k!)) = ^ [l + (27r) d nh(k)} S(k + k'), (C3) 
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where h(k) is given by Eq. f|142|) . Therefore, we obtain the identity 



/(v) 



|e(k,k-v) 



dv = p 1 + (27r) d n/i(k) 



;(k,0) 1 + {27r) d u(k)/3p' 



(C4) 
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